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Abstract. We evaluate the heat conductivity of a dilute pion gas employing the Uehling-Uehlenbeck equa- 
tion and experimental phase-shifts parameterized by means of the SU{2) Inverse Amplitude Method. Our 
results are consistent with previous evaluations. For comparison we also give results for an (unphysical) 
hard sphere gas. 
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1 Transport coefficients 

Transport coefficients in heavy-ion collisions remain of in- 
terest as prospects for direct experimental measurement 
at RHIC improve. In particular, elliptic flow has been pro- 
posed [1] as a tell-tale for viscous effects. There have also 
recently been a number of papers emphasizing viscosity 
(for example [2], [2]) following recent insight from confor- 
mal field theory [3]. Follow-up studies are under prepa- 
ration by several groups from the hadron and from the 
quark-gluon phases in RHIC-theory to see the effect of 
the phase transition on the viscosity. 

Although we have presented a comprehensive study of 
the shear viscosity of a hadronic gas at low temperature 
[5], in this brief report we focuse on another transport 
coefficient, the thermal conductivity, associated to energy 
conservation and flow in the hot gas. 

We employ the same notation and conventions as in 
our earlier publication [5]. In particular we employ the 
SU (2) Inverse Amplitude Method parametrization of the 
pion-pion scattering experimental phase shifts (as well as 
alternative parametrizations to check the sensitivity of the 
calculation presented). Note we also employ the Landau 
convention for the flows as opposed to the Eckart conven- 
tion that has also been widely used . 

To make this paper minimally self-contained, we col- 
lect here a few of the key results and assumptions. 

We are working in a pion gas at temperatures well 
below any phase transition, in the dilute approximation 
near equilibrium. The equilibrium equation of state yields 
an enthalpy per unit volume h = P + p as 
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and the particle (pion) density is given by 
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(z = e^^ ™)/^ being the fugacity). 

We take a constant perturbation around equilibrium 
proportional to the relativistic temperature-pressure gra- 
dient 
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and change the variable to x = [p/m^Y. We perform a 
polynomial expansion oi g{x), and display results for the 
first order only, g{x) = Aq, with a constant. Conver- 
gence is known to be fast. Since pions are not exactly 
massless, no singularities appear that might require a more 
careful treatment and the only reason to improve on 
this variational calculation would be to improve its accu- 
racy beyond the 10% level. This program has been carried 
out at least for the non-relativistic pion gas viscosity in ref. 
[5]. We have also made minimum sensitivity checks based 
on a simple numeric Gram-Schmidt polynomial construc- 
tion beyond Aq but will report them elsewhere. 

From the standard linearized (Boltzmann) Uehling- 
Uehlenbeck transport equation one can derive the follow- 
ing for this perturbation near equilibrium: 
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where w denotes the total energy {uj = E+Ei = E' + E[) , 
A is the inverse normalization constant of the distribution 
functions 
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(ffir = 3 for isospin degeneracy) and the last term is a 
shorthand for the symmctrization 
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Once Aq has been thus computed, the heat conductivity 
follows as 
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To evaluate the integral in I2 in Eq. (|4]) we employ a Mon- 
tecarlo computer program. Without loss of generality we 
can choose the total momentum directed along the OZ 
axis. The independent variables can be taken as P and lo 
(respectively the total momentum and energy in a binary 
collision), p = |p| and p' = |p'| (the incoming and out- 
going pion momenta for one of the two pions, the other 
being fixed by momentum conservation). Finally the out- 
going pion with momentum p' does not need to be in the 
same plane as P and p, and therefore we need an az- 
imuthal angle (f)' for this momentum. The cosines of the 
polar angles associated to p and p' are fixed by the energy 
conservation relation 
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and the associated integrals are immediately performed, 
with 
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In addition, rigid rotations around P parametrized by 4> 
are trivial and lead to a factor of 27r, and global rotations 
of the system (the angles associated to P) are also trivial 
and yield another Att. 

Putting all together, the phase space integration 
weighted with the square scattering amplitude can be ex- 
pressed as 
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Our numerical results for the thermal conductivity are 
displayed in the figures. 

In Fig. [T] we plot the heat conductivity as a function of 
temperature, at several /^'s in the hard-sphere scattering 
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Fig. 1. We show the ^ dependence of the conductivity for the 
hard-sphere gas appoximation 
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Fig. 2. Thermal conductivity k, from SU(2) chiral perturbation 
theory unitarized by means of the Inverse Amplitude Method. 
The pion chemical potential is taken to be ^-n = 0. 



case, that is, when we use a constant scattering amplitude 
based on Weinberg's low energy theorem: 
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For no /.t does the heat conductivity diverge at low tem- 
perature to the reach of our Montecarlo. In the high tem- 
perature limit we expect on dimensional grounds, and nu- 
merically find, K = A ■ T'^ , where the numerical constant 
is close to ^ = 685. 

In Fig. [21 employing lAM phase shifts, that unitarize 
a higher order O(p^) in chiral perturbation theory [9], and 
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largely mediated by resonances (automatically incorpo- 
rated in our Inverse Amplitude Method Phase shifts). 
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Fig. 3. Thermal conductivity times temperature Tn from 
SU(2) chiral perturbation theory unitarized by means of the In- 
verse Amplitude Method. The pion chemical potential is taken 
to be fiTy — 0. We also plot the computation with the phase- 
shifts of TO^ and compare them to D. Davesne's evaluation. 



therefore introduce a dependence with the energy for the 
pion-pion scattering amplitude, we now find that there is a 
minimum value near Tc — 100 MeV and k — > oo when T — > 
0. The high T scaling is close to the dimensional analysis 
even at the moderately large plotted temperatures, 

K (X 3.4T2 . 

In the figure we also show our own Montecarlo-based cal- 
culation of the heat conductivity employing the simple 
resonance saturation parametrization for the isoscalar and 
isovector phase shifts by the Brookhaven group [10], for 
comparison and to give an idea of the sensitivity to the 
parametrization. 

Next in Fig. [3] we show the comparison with the exist- 
ing computation of D. Davesne [TT] of both calculations 

Note also a recent calculation by that employs 
chiral perturbation theory alone (without unitarization). 
This approach features an ever-increasing cross section, 
unphysical behavior that artificially shortens the mean 
free path and therefore lowers the transport coefficients. 
Therefore the validity of the results is limited to low tem- 
peratures. However a direct comparison with this approach 
is difficult since the authors directly include the effect of 
baryon resonances. 

In conclusion, from published calculations and our own 
contribution we see that we have a fair theoretical idea on 
how the thermal conductivity of a pure pion gas should 
behave with temperature. We have further studied its be- 
havior with chemical potential that is important because 
chemical freeze-out is expected to occur before thermal 
freeze out in Heavy-Ion Collisions, the reason being that 
the low energy hadron interactions are elastic scattering. 
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